Neutrino-nucleus reactions in the delta 
resonance region 



B. Szczerbinska a ' c T. Sato b K. Kubodera c T. -S. H. Lee d 

a Dakota State University, Madison, South Dakota 57020, USA 

b Department of Physics, Osaka University, Toyonaka, Osaka 560-0043, Japan 

c Department of Physics and Astronomy, University of South Carolina, Columbia, 

South Carolina 29208, USA 

d Physics Division, Argonne National Laboratory, Argonne, Illinois 604-39, USA 



Abstract 

Reliable estimates of neutrino-nucleus reactions in the resonance-excitation region 
play an important role in many of the on-going and planned neutrino oscillation ex- 
periments. We study here neutrino-nucleus reactions in the delta-particle excitation 
region with the use of neutrino pion-production amplitudes calculated in a for- 
malism in which the resonance contributions and the background amplitudes are 
treated on the same footing. Our approach leads to the neutrino-nucleus reaction 
cross sections that are significantly different from those obtained in the conventional 
approach wherein only the pure resonance amplitudes are taken into account. To 
assess the reliability of our formalism, we calculate the electron-nucleus scattering 
cross sections in the same theoretical framework; the calculated cross sections agree 
reasonably well with the existing data. 
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1 Introduction 

It is well recognized that the precise knowledge of neutrino-nucleus reaction 
cross sections is of importance in analyzing neutrino oscillation experiments; 
for recent reports, see e.g. Refs. [1,2,3,4,5]. In particular, neutrino-nucleus 
reactions at incident neutrino energies around 1 GeV play a prominent role 
in many cases including the experiments at K2K [2]. To obtain estimates of 
the relevant cross sections, one must at present rely on theory, and much 
theoretical effort has been invested to provide these estimates [6,7,8,9,10,11]. 
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In an attempt to make a quantitative estimation of neutrino-nucleus reaction 
cross sections, it is useful to study simultaneously the related electron-nucleus 
reactions within the same general theoretical framework, and this strategy 
has been pursued by many authors. In electron-nucleus scattering in the GeV 
region, quasi-elastic scattering and pion-production processes are known to be 
the main reaction mechanisms, and similar features are expected to manifest 
themselves also in the neutrino-nucleus reactions in the GeV region. 

For quasi-elastic scattering, the relevant transition operators are essentially 
known, so the main theoretical issue is how to incorporate various nuclear 
effects for the initial and final states. The early works were based on the 
Fermi gas model [12,13], but recent investigations incorporate the nuclear 
correlation effects in the initial state with the use of the spectral function and 
take account of the final-state interactions on the outgoing nucleon [7]. As 
regards the pion-production process, in addition to these nuclear effects, the 
structure of the transition operators responsible for pion production needs to 
be carefully studied. These operators can in principle involve more than one 
nucleon, but it is in general expected that pion production on a single nucleon 
should give a dominant contribution. Neutrino-induced pion production on the 
nucleon in the resonance region has been studied so far mostly with the use of 
pure resonance excitation amplitudes. In some studies these amplitudes were 
evaluated in the quark model, see, e.g., Ref. [14]. In recent studies by Paschos 
and his collaborators [15,16], the resonance excitation amplitudes due to the 
vector current were directly related to the empirically known electro-excitation 
amplitudes, while those due to the axial-vector current were constrained by 
invoking PCAC. 

Meanwhile, it is to be noted that pion production can take place not only 
through resonance excitations but also via non-resonant processes. Two of 
the present authors [17,18] have recently developed a dynamical model for 
describing photo- and electro-production of pions off the nucleon around the 
A-resonance region, with the view to systematically incorporating both the 
resonance and non-resonance contributions. Hereafter we refer to this approach 
as the SL-model (the Sato-Lee model). The development of the SL-model was 
motivated by recent extensive experimental studies of electron- and photon- 
induced meson-production reactions on the nucleon in the resonance region. 
The main objective of these experiments is to study the non-perturbative 
features of QCD by testing the resonance properties as predicted by QCD- 
inspired models and/or lattice simulations. The SL-model was subsequently 
extended to weak-interaction processes [19,20], and it was shown that this 
model gives a successful description of neutrino-induced pion production in 
the A-resonance region. 

As explained in more detail later, the SL-model starts from the non-resonant 
meson-baryon interaction and the resonance interaction, and the unitary am- 
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plitudes are obtained from the scattering equation. It leads to fairly consis- 
tent descriptions of all the available data for the electroweak reactions in the 
A-resonance region. It has been shown that treating the resonance and non- 
resonance amplitudes on the same footing can have significant observable con- 
sequences. In particular, the inclusion of the pion cloud effects as considered in 
SL can resolve a long-standing puzzle that the N-A magnetic dipole transition 
form factor Gm predicted by the quark model is smaller than the empirical 
value by as much as ~40%. Furthermore, the electric E2(G S ) and Coulomb 
C2(G C ) form factors for the N-A transition in electron scattering calculated 
in the SL-model show pronounced momentum dependences due to the pion 
cloud effects, which suggests non-negligible deformation effects in the iV-A 
transition. Regarding the neutrino reactions, a serious problem that has been 
known for quite some time is that the axial-vector N-A transition strength 
calculated in the constituent quark model [21] is lower than the empirical 
value [22] by about 35 %. It is noteworthy that the dynamical pion cloud 
effects included in the SL-model [19] can naturally remove this discrepancy. 

In view of these successes, it seems worthwhile to study neutrino-nucleus re- 
actions in the resonance region with the use of the SL-model amplitudes for 
neutrino-induced pion production on the nucleon. We describe here our first 
attempt at such a study and present the cross sections, the energy spectrum 
of the final lepton (for charged-current reactions), and the lepton- momentum 
transfer distribution. Our work is basically of exploratory nature and, as far 
as the nuclear effects are concerned, we only consider those that can be taken 
into account with the use of a modified Fermi gas model wherein nuclear cor- 
relations are approximately subsumed into the spectral function [7]. Despite 
these limitations, our investigation is hoped to be informative as the first cal- 
culation of neutrino-nucleus reactions in the A-resonance region based on the 
electroweak pion-production amplitudes calculated in SL [17,18], whose va- 
lidity has been extensively tested by the Jlab data [23,24]. It is understood 
that, as the experimental precision improves, more detailed calculations will 
be called for that incorporate higher order effects. In particular, the final-state 
interaction (FSI) must be treated properly. As discussed in Ref. [7] and many 
earlier works on inclusive electron scattering, FSI re-distributes the inclusive 
cross sections and, for the incident electron energy around 1 GeV, FSI can 
reduce the strength at the quasi-free peak by about 10 %. We remark that a 
detailed study [25] indicates that the FSI effects for inclusive reactions, prop- 
erly treated, can be rather different from those for exclusive processes [26]. 
Meanwhile, in the pion production region, we need to take into account pion 
absorption and medium effects on A propagation. To this end, one may prof- 
itably use the information obtained in the well-developed A-hole model [27]; 
such a study has been made in Ref. [28] within the framework of the dynamical 
transport approach. Our present calculation, however, falls short of consider- 
ing FSI. Since the importance of these FSI effects grows rather fast with the 
increasing target mass number (this is particularly true for pion absorption), 
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we limit ourselves here to nuclear targets of low mass numbers and concentrate 
on the A=12 target. Although heavier nuclei such as 56 Fe are important in 
some neutrino-oscillation experiments [5], we can deal with these cases only 
after FSI is incorporated into our formalism. 



2 Sato-Lee (SL) Model 



As the SL-model has been fully described in Refs. [17,18,19,20], we give here 
only a brief explanation of the model, using as an example the case of pion 
photoproduction. The effective Hamiltonian H e ff in the SL-model for this 
process is given by 



Heff — Hq + V wN + f 77r + T wN ^ A + r 7A r^ A , (1) 

where H Q is the free Hamiltonian; i^jv and t> 7vr represent the non-resonant pion- 
nucleon and pion photoproduction interactions, respectively, while T^^a and 
r 7 jv<-^A are responsible for the creation and annihilation of a bare A-resonance. 
By solving the Lippmann-Schwinger equation based on the above effective 
Hamiltonian, we obtain the amplitude for pion production on a nucleon as 



where E is the total energy of the pion and nucleon in the center-of-mass 
system. The first term t 77r is the non-resonant amplitude, which arises from 
the vertices v n N and v 7 tv alone, while the second term represents the resonant 
amplitude involving the dressed vertex f . We note that the bare resonance 
vertex T is renormalized into f by the non-resonant meson cloud effects arising 
from rescattering as 



r 7 AT_ > A(-E) — r 7 Ar_»A + ^wN^AGotyn(E), (3) 

and that the renormalized resonance vertex f should exhibit a significant 
deviation from the bare vertex T because of the meson cloud effects. In con- 
ventional analyses, however, one disregards the difference between the bare 
and renormalized resonance vertices and, assuming the Breit-Wigner form for 
(what in our approach is identified as) the renormalized vertex, tries to extract 
the parameters characterizing that form by fitting to the data. However, as 
discussed in detail in Refs. [17,18,19,20] and as briefly mentioned in the intro- 
duction, the resonance properties deduced from this simplified treatment tend 
to exhibit significant discrepancies with the theoretical predictions, indicating 
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the importance of considering the resonant and non-resonant contributions 
simultaneously. 

We expect that the unified treatment of the resonant and non-resonant con- 
tributions should have significant consequences in the neutrino-nucleon and 
neutrino- nucleus reactions as well. To illustrate this point, we give in Fig. 1 
the electron energy spectrum da/dE e for the u e N — > e^irN reaction calculated 
in the SL-model. Since the isospin of the 7rN system for the u e p — > e~n + p reac- 
tion is 3/2, one may naively expect that the A-resonance amplitude dominates 
the cross section. However, the energy spectra obtained in our SL-model cal- 
culation (shown in the solid lines) are markedly different from those obtained 
from the dressed resonant amplitudes (shown in the long-dash lines). For com- 
parison, the results obtained in the Lalakulich-Paschos (LP) model [15] are also 
shown in the short dashed lines. If in our approach we drop the contribution 
of the non-resonant amplitudes (retaining only the resonance contributions), 
the results turn out to be very similar to those obtained in the LP-model. 




Fig. 1. Electron energy spectrum for v e p — > e~ir + p (left panel), u e n — > e~7r + n 
(middle panel), and u e n — > e~ir°p (right panel). For the explanation of the three 
curves in each panel, see the text. 

Since neutrino-nucleus reactions obviously involve both up and un reactions, 
it is also informative to examine how the non-resonant contributions can af- 
fect the relative importance of the up and un contributions. If the 7 = 3/2 
resonance amplitude dominates, the cross section on the neutron should be 
1/3 of that for the proton. Fig. 2 shows da/dQ 2 for the proton target (solid 
line) and 3 x da/dQ 2 for the neutron target (long-dashed line) for E v = 1 
GeV. The curves should agree with each other if the delta mechanism domi- 
nates. However, the neutron cross section is about 20% larger than the value 
expected from A dominance. The short-dashed line gives da/dQ 2 for the pro- 
ton target obtained with the use of the dressed resonance amplitude alone. 
We note that this curve overlaps rather well with the solid line correspond- 
ing to the SL-model results, a feature to be contrasted with the behavior of 
da/dEi shown in Fig. 1. We learn from this that, even for the same reaction, 
some observables are more sensitive than others to the difference between the 
SL-model results and those of the resonance-contribution-only approach. 
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Fig. 2. Differential cross sections da/dQ 2 for v e N — > e~irN at E v = 1 GeV, where 
Q 2 = — g 2 and q = p e — p u is the lepton momentum transfer. Solid line - proton tar- 
get (SL-model calculation); short-dashed line - proton target (resonance amplitude 
only); long-dashed line - neutron target case multiplied with a factor of 3 (SL-model 
calculation). The final irN state is n + p for the proton target, whereas both and 
ir + n can contribute in the neutron target case. 

3 Neutrino-nucleus reaction 



We consider the charged-current (CC) neutrino-nucleus reaction 



v e (p v ) + \i(Pi) >^ i{pt) + \f(P f ) > , (4) 

where \i > represents a target nucleus of mass A, \f > stands for a final 
hadronic state, and £ is a lepton flavor (£ — e, //,r); the relevant momenta 
are indicated in the parentheses. The cross section for this process is given in 
terms of the lepton tensor and the hadron tensor W^ u as 



da _ Pi G% cos 91 (5) 



dEtdtti Pu 8tt 2 
The lepton tensor is expressed as 

IT = 2 [pfp: + rftf - gr{p v ■ Vi - m 2 ) + ie^p„, aPl , p ] , (6) 

where me is the mass of the final lepton. The hadron tensor is related to the 
matrix elements of the hadronic weak current J M as 



w = E(2vr)V5 4 (p / + Pl - Pi - Vv ) < f\j»\i>< f\r\i >*, 

f 
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where V is the quantization volume. In the Fermi-gas model, W^ u can be 
related to the single-nucleon transition amplitudes, and the relation for the 
case of quasi-elastic scattering is well known [12] . For a single-pion production 
process, W^ u is given as 

W vv = J dptdkdp6(p F -\p\)6(\fi\ - p F ) 5\p + q - p' - k) 

1 m\ 3 

X (2vr) 3 E N {p)E N {p')2E^k)Anp F 

AT, 

x E ~Y~ <K l N(p', s N ,,t N ,) \f\N(p, s N , t N )> 

x <7r i iV(p', Sjv ,,t JV ,)|j1^(p, Sj v,t J v)>* (7) 

Here p (pf) is the four-momentum of the initial (final) nucleon, k is the four- 
momentum of the pion, q = p v — pi, and p F is the Fermi momentum; N tN = 
Z/N are the proton and neutron numbers in the target nucleus; sjy and 
are the spin and isospin of the struck nucleon, while % is the isospin index 
of the pion. The matrix element of the nucleon current for pion production, 
<7r' l N(p', SN',tN')\j^\N(p, SN,t N ) >, is calculated using the SL model. Since 
SL gives the pion-production amplitude in the pion-nucleon center-of-mass 
frame (7riV-cm frame, for short), we transform it into the amplitude in the 
rest frame of the target nucleus (LAB frame) according to 



s N ,s N /,t N ,t 

xN tN Jdnj(\pf\- PF ) 



m N 



, Anp%J y 1 u E N (p) 

S N ,S N l,t N ,l ff 

\kc\rriN 
32tt 2 W 

x < 7rW(p', s N .) \f \N(p, s N , t N ) > 

x A^'<n*N(p', s N ,)\f'\N(p, s N , t N ) >l N ^ cm , (8) 



where W is the invariant mass of the pion and nucleon given 

The Lorentz transformation matrix transforms vectors in the 7riV-cm 

frame to those in the LAB frame. In the 7rA^-cm frame, p'+k = (W, 0) , whereas 



in the LAB frame we identify p' + k = p + q = (yp 2 + mjy — B + u,p + q); 
thus the nuclear binding correction is taken into account with the use of 



= (0O 2 + m 2 N — B,p). We note that the Pauli blocking factor, 9(\pf \ — p F ), 

— * — * 

is dependent on the pion momentum k through pi = p + q — k, and hence 
the consideration of the Pauli blocking effect requires the knowledge of the 
pion-production amplitude. We come back to this point later. 

We take into account the nuclear correlation effects in the initial state by 
using the spectral function P(p ) E) obtained in Ref. [29]. This is achieved by 
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the following replacements in Eq. (8) 



^3- J dp9(PF-\p\) - / dpdEP{p, E) (9) 
P° = \jv> + m 2 N -B^m N -E. (10) 

We note that P(p,E) is normalized as JdpdEP(p, E) = 1. Although the use 
of the spectral function implies that the separation of occupied and empty 
nucleon orbits based on the Fermi momentum pp is no longer strictly valid, 
we choose to retain the factor 9(\pt\ — Pf) in Eq.(8) to approximately take 
account of Pauli blocking for the final nucleon. 



4 Results and Discussion 



Using the formalism explained in the previous section, we calculate neutrino- 
nucleus reaction cross sections for a representative case of the z/- 12 C scatter- 
ing. For the incident neutrino energy we take E v — 1 GeV, a value lying 
in the energy region of current importance for many neutrino oscillation ex- 
periments. We first discuss the pion-production cross sections, which are our 
main results, and subsequently we consider the combined contributions of the 
pion-production and quasi-elastic processes. 

The differential cross sections for the v X2 Q,-^ e~-nX reaction, normalized with 
the target mass number (A = 12), are shown in Fig. 3, for the lepton scat- 
tering angle 6 = 10° and 30°, as a function of the invariant mass W = 
\j —Q 2 + m 2 N + 2urniN- The dashed curve is the cross section obtained sim- 
ply by taking the average of the incoherent contributions of the free protons 
and neutrons, while the Fermi-gas model results are shown by the dash-double- 
dotted curves. As expected, the inclusion of the nucleon Fermi motion widens 
the resonance width compared with the free nucleon case. The dash-dotted 
curves, corresponding to the case that includes the Pauli blocking effect for 
the final nucleon, indicate that the blocking effect reduces the forward cross 
section by about 20%. As mentioned in connection to Eq. (8), the inclusion 
of the Pauli blocking effect for the 7r-production process requires the knowl- 
edge of the pion-production amplitude. This implies that this effect cannot be 
evaluated by taking (as often done in the literature) the incoherent sum of the 
free-nucleon pion-production strength over the Fermi sea. (The role of Pauli 
blocking, however, diminishes for larger angles, where the momentum transfer 
becomes larger than the Fermi momentum.) The solid curves show the results 
of our full calculation that includes the spectral function taken from Ref. [29]; 
it is seen that the nuclear correlation effects further broaden the peak width 
and reduce the peak height by about 20%. 
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Fig. 3. Differential cross sections for vg C— > e"7rX at 6 = 10° and 30°. Dashed line 
- free-nucleon; dashed-double-dotted line - Fermi-gas model; dashed-dotted line - 
Fermi-gas model with Pauli blocking; solid line - full calculation. 




Fig. 4. Differential cross sections for i/ e 12 C — > e~7rJf at = 10° and 30°. The dashed 
line represents the case where the pion production amplitude contains the dressed 
A contribution alone, while the solid line shows the results of our full calculation. 

As discussed, the non-resonant mechanism contained in the SL-model plays 
a more pronounced role for the neutron than for the proton. To illustrate 
how this feature affects the neutrino-nucleus reaction, we give in Fig. 4 the 
differential cross sections for z/g 2 C-^ e~nX calculated with and without the 
non-resonant contributions; the solid curve is the result of the full calculation, 
while the dashed curve presents the case where only the contribution of the 
dressed resonance amplitude is considered. As expected, the resonance-only 
approach underestimates the cross sections by about 20% even in the reso- 
nance region. The full calculation is found to give more strength for lower 
values of W than the resonance-only case. 
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Fig. 5. Left panel - Differential cross section for the i/g C— > e~X reaction at E v = 1 
GeV and 9 = 30°. Right panel - Differential cross section for the e~ + 12 C — > e~V 
reaction at i£ e = 1.1 GeV with e = 37.5°. The experimental data points are from 
Ref. [30]. 



To the contribution of the pion-production process we now add the contri- 
bution of the quasi-free nucleon knockout process. The latter is calculated 
using again the modified Fermi-gas model that incorporates the spectral func- 
tion [7,29]. At the incident energy under consideration, the sum of these two 
contributions is expected to give the bulk of the inclusive reaction cross sec- 
tion. The differential cross section for z/g 2 C— > e~X at 9 = 30° is shown in the 
left panel in Fig. 5. The bump at the lower energy is due to quasi-free nucleon 
knockout, while the higher energy bump is due to A- resonance excitation. To 
examine the validity of our present approach, we apply the same calculational 
framework to the e~+ 12 C — > e~X reaction (with the weak current replaced by 
the electromagnetic current), and compare the results with the experimental 
data. Fig. 5 shows this comparison. It is seen that the general trend of the data 
is reproduced reasonably well; in particular, the magnitude of the cross section 
in the A-resonance region is well reproduced. We remark that a calculation by 
Benhar et al. [7](BFNSS) underestimates the height of the A peak. According 
to Ref. [31], this is perhaps mainly due to the fact that, around the A region, 
the neutron structure functions of Ref. [32] used in BFNSS are significantly 
weaker than those extracted from an analysis of inclusive electron scattering 
on the deuteron. It is noteworthy that the SL model, in addition to providing 
a satisfactory description of the proton structure functions [20] , gives neutron 
structure functions that are in good agreement with those deduced in Ref. [31]; 
see Fig. 6. This feature explains the difference in the A peak height between 
our results and those of BFNSS. 



We note, however, that our calculation gives a dip structure that is somewhat 
too deep, a feature that seems to indicate that we need to go beyond 'impulse' 
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Fig. 6. The neutron structure function W 2 n at E e = 2A45GeV and e = 20°. The 
solid line represents the results of the SL-model and the shaded area represents the 
W£ from the analysis of Ref. [31]. 

approximation and/or employ more elaborate treatments of nuclear correla- 
tion effects; see also Ref. [9] . It is also to be noted that in the higher W region 
our model, which only includes the A-resonance, is likely to underestimate the 
transition strength. 

For some purposes it seems useful to present our results in the form of Q 2 - 
distribution [Q 2 = — (p v — pi) 2 ] or ^-distribution. We get da/dQ 2 and da/dEi 
using the formulas, 



Fig. 7 gives the Q 2 and spectra for the v l2 Q reaction at E v = 1 GeV; the 
left (right) panel corresponds to the CC (NC) reaction. The total contribution 
(solid line) consists of the quasi-free contribution (dash-dotted line) and the 
pion production contribution (short-dashed line). We note that the pion pro- 
duction contribution is reduced by the introduction of the spectral function. 
Finally, the muon energy distribution da/dE^ in the v X2 G— > [TX reaction is 
shown in Fig. 8. 



5 Summary 

We have studied the neutrino-nucleus and electron-nucleus reactions in the 
A-resonance region with the use of the Sato-Lee (SL) model, which allows us 
to treat the resonant and non-resonant contributions in a unified manner in 
deriving the amplitudes for pion electroweak production on a nucleon. The 
validity of SL has been extensively tested for electromagnetic observables in- 
volving single-nucleon targets, and we can use the available electron-nucleus 
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Fig. 7. The differential cross section da/dQ 2 for ujfC^ and u 12 C^ vX at 

E v =1 GeV. The sum of the quasi- free and pion-production contributions is shown 
by the solid line (full calculation), and by the long-dashed line (Fermi-gas model). 
The individual contribution of the quasi-free process is shown by the dash-dotted 
line, and that of the pion-production process by the short-dashed line. 
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Fig. 8. Muon energy distribution da/dE^ for the vj^C— > pTX reaction. The con- 
tributions of the quasi-free and pion-production reactions are shown in the short- 
-dashed and dash-dotted lines, respectively, while their sum is given by the solid 
line. The long-dashed line represents the results obtained in the Fermi-gas model. 

scattering data to assess the reliability of the application of SL to neutrino- 
nucleus reactions. As for the nuclear correlation effects, we have considered 
here only those effects which can be considered to be subsumed in the spec- 
tral function. Despite this rather limited treatment of the nuclear effects, our 
calculation based on SL gives reasonably good descriptions of the relevant elec- 
tron scattering data; the peak structure in the cross section in the resonance 
region is well reproduced by our calculation. It is reasonable to expect that 
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our calculation of the neutrino-nucleus reaction cross sections based on the 
SL model enjoys the same level of success. It seems worthwhile to further de- 
velop SL studies of neutrino- and electron-nucleus reactions by elaborating the 
treatment of the nuclear effects (including medium effects on the A-resonance 
itself). 
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